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Abstract 

Ozone is an important radiative trace gas in the Earth’s atmosphere and has also 
been detected on Venus and Mars. The presence of ozone can significantly influence 
the thermal structure of an atmosphere due to absorption of stellar UV radiation, and 
by this e.g. cloud formation. Photochemical studies suggest that ozone can form in 
carbon dioxide-rich atmospheres. Therefore, we investigate the effect of ozone on the 
temperature structure of simulated early Martian atmospheres. 

With a ID radiative-convective model, we calculate temperature-pressure profiles 
for a 1 bar carbon dioxide atmosphere containing various amounts of ozone. These 
ozone profiles are fixed, parameterized profiles. We vary the location of the ozone 
layer maximum and the concentration at this maximum. The maximum is placed at 
different pressure levels in the upper and middle atmosphere (1-10 mbar). 

Results suggest that the impact of ozone on surface temperatures is relatively small. 
However, the planetary albedo significantly decreases at large ozone concentrations. 
Throughout the middle and upper atmospheres, temperatures increase upon introduc¬ 
ing ozone due to strong UV absorption. This heating of the middle atmosphere strongly 
reduces the zone of carbon dioxide condensation, hence the potential formation of car¬ 
bon dioxide clouds. For high ozone concentrations, the formation of carbon dioxide 
clouds is inhibited in the entire atmosphere. In addition, due to the heating of the 
middle atmosphere, the cold trap is located at increasingly higher pressures when in¬ 
creasing ozone. This leads to wetter stratospheres hence might increase water loss 
rates on early Mars. However, increased stratospheric H 2 O would lead to more HO^, 
which could efficiently destroy ozone by catalytic cycles, essentially self-limiting the 
increase of ozone. This result emphasizes the need for consistent climate-chemistry 
calculations to assess the feedback between temperature structure, water content and 
ozone chemistry. Furthermore, convection is inhibited at high ozone amounts, leading 
to a stably stratified atmosphere. 


* Corresponding author: email vonparis@obs.u-bordeauxl.fr, tel. -1-33 (0)557 77 6131 


Preprint submitted to Elsevier 


May 29, 2015 



Keywords: 

early Mars; habitability, atmospheres 


1. Introduction 


The climate and, consequently, the habitability of early Mars are a long-standing 
question that has not yet been answered conclusively. From an atmospheric modeling 
point of view, short episodes of a warm, wet climate, followed by long periods of arid, 
cold climates, seem to be favored for the Noachian period about 3.8 billion years ago 
(e.g., [Segura et al.|2008|[Wordsworth et al.|2013|[Halevy and Head|20T4l l. 

An important input for atmospheric modeling studies is the atmospheric compo¬ 
sition, together with the surface pressure. Most previous studies (e.g., Kasting[ri991 


Wordsworth et al.|2013| l assumed mixed carbon dioxide (C 02 )-water vapor (H 2 O) at¬ 
mospheres. In these studies, the partial pressure of carbon dioxide was of the order of 
bars, consistent with outgassing model studies (e.g., Phillips et al.pOOT] [Grott et ah] 
201 Ij l and in-situ analyses (e.g.. Manga et al.|2012 [Kite et al.|2014| l. In ID models, the 
atmospheric H 2 O content is usually controlled by a fixed relative humidity to simulate 
a hydrological cycle. 

However, other atmospheric (trace) gases such as molecular nitrogen (N 2 ), molec¬ 
ular hydrogen (H 2 ), methane (CH 4 ) or sulphur dioxide (SO 2 ) might have been present 
in the early Mars atmosphere. The initial N 2 inventory is thought to be relatively large 
(e.g., [McKay and Stoker|[T989|l and can provide modest surface warming (von Paris 


and 


et al. 2013b| l. Methane has also been suggested to warm the atmosphere (e.g., [Postawko 
Cuhn[[1986 Ramirez et al.|2014)l, but the effect on surface temperature is small. 


Ramirez et al. ( 2014| l find that H 2 -induced warming could raise surface temperatures 
above freezing, assuming H 2 was a major atmospheric constituent (around 20 % vol¬ 
ume mixing ratio). The effect of SO 2 has been investigated by a number of studies 


(e.g., PostawkoandKuhn|l9^ Yung et al.|1997 Halevy et al.[[2007 Johnson et al. 


2008 Tian et al.|2010 Mischna et al.|2013] l, but results suggest that it most likely does 

not contribute strongly to warming the surface. 

One proposed solution to warming early Mars has been the formation of carbon 
dioxide clouds (e.g., Pierrehumbert and Erlick|1998 [Forget and Pierrehumbert|1997j l. 
Early ID modeling studies suggested that the needed warming however strongly de¬ 


pends on the assumed cloud cover, which would have to be nearly 100% (e.g., Mischna 
e^jJj 2000 ) 1 . Following time-dependent ID modeling studies by Colaprete and Toon 
( 2003 [ ) found this to be unrealistic, a conclusion also supported by recent 3D studies 
(e.g.. Forget et al.|2013[ [Wordsworth et al.|2013| l. In addition, [Kitzmann et al.[p013[ l 
suggested that previous radiative transfer algorithms probably strongly overestimated 
the warming associated with carbon dioxide clouds. 

The formation of carbon dioxide clouds depends on the temperature profile in the 
middle atmosphere. As shown by, e.g., Yung et al. \\991\ or Ramirez et al. P014| l, UV 
absorption by SO 2 or near-IR absorption by CH 4 can effectively inhibit carbon dioxide 
condensation, or at least reduce the cloud formation region. 

So far, ozone (O 3 ) has not been considered in the context of early Mars. Like 
SO 2 , ozone has strong UV and visible absorption bands. On Earth, these are responsi- 
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ble for the pronounced stratospheric temperature maximum. Ozone has been detected 
in the C 02 -rich atmospheres of both Venus (e.g., Montmessin et al.||2011|l and Mars 


(e.g., [Lebonnois et al.||2006ll. Furthermore, numerous photochemical modeling stud- 


ies suggests that ozone can be formed from abiotically produced oxygen (e.g., Sel- 

sis et al .|2002 

Segura et al.|2007 Domagal-Goldman and Meadows|2010 

|Domagal- 

Goldman et a 

2014|l. In addition, Wordsworth and Pierrehumbert 

m\A 

1 propose a 


(non-chemical) abiotic oxygen formation process based on hydrogen escape and water 
cold-trapping. Therefore, it is reasonable to assume that on early Mars, ozone could 
have been also, to a certain extent, present in the atmosphere. Hence, in this study, we 
explore the effect of ozone on the temperature structure of early Mars 

The paper is structured as follows: Section [^presents the atmospheric model and 
the simulations performed. Results are shown and discussed in Sect, [^and conclusions 
in Sect. |4] 


2. Computational details 

2.7. Model description 

We use a ID, steady-state, cloud-free radiative-convective atmosphere model to 
calculate globally, diurnally averaged temperature and H 2 O prohles. The model is 
originally based on Kasting et al. ( 1984a| l and Kasting et al. (1984b 1 . Further code 
developments are described in, e.g., |Kasting| (|1988 1 , Mischna et al. (2000|l, |von Paris 

aLjjSoToF 


|et al.| ( |2008| l or |von Paris et al.| 

The model atmospheres are divided into 52 levels. Temperature prohles in the up¬ 
per atmosphere are obtained by solving the radiative transfer equation (38 bands for 
incoming stellar radiation, 25 bands for planetary and atmospheric thermal radiation). 
For incoming stellar radiation (0.2-4.5 pm), Rayleigh scattering (by N 2 , H 2 O, CO 2 and 
CH 4 ) and molecular absorption (by H 2 O, CO 2 , CH 4 and O 3 ) contribute to the opac¬ 
ity. To allow for scattering, the angular integration of the radiative transfer equation is 
performed using a 2-stream code ( |Toon et al.|1989| l. For thermal radiation (1-500/rm), 
molecular absorption (by H 2 O, CO 2 , CH 4 and O 3 , see below. Sect. 2.2 1 and contin¬ 
uum absorption (by N 2 , H 2 O and CO 2 ) are considered. N 2 continuum data is taken 


from Borysow and Frommhold ( 1986[ ) and |Lafferty et al. ( |1996| l, whereas H 2 O self 
and foreign continua as well as the CO 2 foreign continuum are incorporated following 
Clough et al. ( 1989|l. The CO 2 collision-induced self continuum is described follow¬ 
ing Kasting et al. ( 1984b| ). In the lower atmosphere, the model performs convective 
adjustment such that temperature prohles follow the wet adiabat. The formulation of 


the adiabatic lapse rate takes into account the condensation of H 2 O or CO 2 (Kasting 
]V9]^ The super-saturation ratio for the onset of CO 2 condensation is unity (but see 


Glandorf et al.| ( 2002] l for a different estimate for early Mars). 

H 2 O prohles are re-calculated for each model time step, according to local tem¬ 
perature and using a hxed relative humidity (RH) prohle. The concentrations of the 
other species are hxed at the start of the calculations and only adjusted according to 
condensation of water and carbon dioxide at the surface ( |von Paris et al.|2013a| l. 

For more details, we refer to|von Paris et al.|(|2008|l and references therein. 
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Table 1: Spectral intervals for the new IR radiative transfer scheme and species considered 

Interval range [cm“'] contributing species 


1 

7,470 - 

10,000 

C02,H20,CH4 

2 

6,970 

- 7,470 

C02,H20,CH4 

3 

6,000 

- 6,970 

C02,H20,CH4 

4 

5,350 

- 6,000 

C02,H20,CH4 

5 

4,600 

- 5,350 

C02,H20,CH4 

6 

4,100 

- 4,600 

C02,H20,CH4 

7 

3,750 

-4,100 

C02,H20,CH4 

8 

3,390 

- 3,750 

C02,H20,CH4 

9 

3,050 

- 3,390 

C02,H20,CH4 

10 

2,750 

- 3,050 

C02,H20,CH4 

11 

2,400 

- 2,750 

C02,H20,CH4 

12 

2,250 

- 2,400 

C02,H20,CH4 

13 

2,150 

- 2,250 

C02,H20,CH4,03 

14 

2,000 

-2,150 

C02,H20,CH4,03 

15 

1,850 

- 2,000 

C02,H20,CH4,03 

16 

1,400 

- 1,850 

C02,H20,CH4,03 

17 

1,100 

- 1,400 

C02,H20,CH4,03 

18 

1,000 

- 1,100 

C02,H20,CH4,03 

19 

905 - 

1,000 

C02,H20,CH4,03 

20 

820 

-905 

C02,H20,CH4,03 

21 

730 

-820 

C02,H20,03 

22 

600 

-730 

C02,H20,03 

23 

525 

-600 

C02,H20,CH4,03 

24 

460 

-525 

C02,H20,CH4 

25 

20 - 

460 

H20,CH4,03 


2.2. Model improvements 


The original IR radiative transfer scheme MRAC (as used in, e.g., von Paris et al. 


2008| von Paris et al.|20i0 von Paris et al.|2013b|l only considered gaseous absorption 


by H 2 O and CO 2 for the calculation of opacities. For this work, a new, updated MRAC 
version has been created that also considers CH 4 and O 3 for gaseous absorption. 

Table [T] shows the spectral intervals as well as the species included in these inter¬ 
vals. 

Line positions and line strengths for both CH 4 and O 3 have been taken from the 
Hitran 2008 database (Rothman et al.|[2()09|l. Line strengths are converted from the 


reference temperature of 296 K to the desired temperature following Norton and Rins- 


|land| ( [T99T| ). Cross sections are then obtained with the MIRART-Squirrl line-by-line 
radiative transfer code ( Schreier and Schimpf|2001 Schreier and Bottger|2003 1 using 
10^ spectral points per band. The line shape is taken to be a Voigt line profile with a 
cutoff of 10 cm '. The foreign component of the Lorentz broadening is taken directly 
from the Hitran database, i.e. for air. The self component of the Lorentz broadening 
was calculated assuming a volume mixing ratio (vmr) of 10 ® for both CH 4 and O 3 , 
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Table 2: Interpolation and range of the grid variables in the new IR radiative transfer scheme 


Quantity 

Temperature 

Pressure 

H 2 O concentration 
carbon dioxide concentration 
CH 4 concentration 
O 3 concentration 


Range Interpolation 
linear in T 
linear in log(p) 
linear in vmr 
linear in vmr 
linear in vmr 
linear in vmr 


Comments 

AT=50K 

Alog(p)= 0.1 

1 -4 points per order of magnitude 
1 -4 points per order of magnitude 
1 -4 points per order of magnitude 
1 -4 points per order of magnitude 


100-400 K 
10-5-1.5 bar 
lO-"^-! 
10-5-1 
io-^-io-2 
io-^-io-2 


3.55 10-^ for CO 2 and 10-^ for water. Numerical tests showed that except for high 
water concentrations (> 10 -^), the effect of the self component on the cross sections 
was negligible. 

The chosen line cutoff in this work is relatively short. Tests with a line-by-line 
radiative transfer code with line cut-offs at 25 cm-^ have shown that the overall flux 
change rarely exceeds 5 % in a given spectral band. For the total integrated outgo¬ 
ing thermal flux, we And a difference of about 1 %, equal to a radiative forcing of a 
few W m-^. This could lead to a small change in surface temperature of a few K (see 
also, for example, Wordsworth et al.|2010 l. However, this is not expected to qualita¬ 
tively alter our conclusions. Furthermore, as pointed out by, e.g., |Halevy et al.| (|2009|l, 
Wordsworth et al. (20101 or Mischna et al. ( 2012| l, assuming a Voigt line profile is 
probably not a well-justified choice in dense CO 2 atmospheres where the far wings of 
lines are substantially sub-Lorentzian. Sensitivity tests by Wordsworth et al. (20T0| or 
von Paris et al. ( 2010| l showed that the influence on surface temperature can be large (of 
the order of 10 K). Still, this is again not likely to influence the conclusions regarding 
the impact of ozone on the thermal strTicture in the upper and mid atmosphere. 

The interpolation variables are temperature, log(pressure) and the individual con¬ 
centrations c, of the absorbing species. The concentration grid has four dimensions, 
one for H 2 O, CO 2 , CH 4 and O 3 , respectively. Table summarizes the range of the 
interpolation variables and the type of interpolation. 

Figure illustrates our approach to calculate the needed k distributions. In every 
spectral interval, for each temperature-pressure point (T/p), cross sections cr, have been 
calculated for the contributing species i (see Table [TJ. These are then combined to 
obtain effective cross sections o-^ff for the different gas mixtures according to assumed 
molecular concentrations c, from Table 


CTefF 




( 1 ) 


From the effective cross sections crgff, k distributions have been obtained that are 
then used in the model for the radiative transfer calculations with a cotTelated-k ap¬ 
proach. 


2.3. Model verification 

We tested the new model against results for early Mars from 
As an example, Figure|^shows the temperature profiles and differences for a 1 bar car- 


von 


Paris et al. ( 2013b| l. 
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Figure 1: Flowchail for calculation of k distributions. 


bon dioxide case (without N 2 ). It is clearly seen that both models compare relatively 
well with each other. The small discrepancies are due to a slightly different approach to 
calculate the line broadening parameters for the k distributions at high pressures, com¬ 
pared to earlier work (e.g., von Paris et al.|[2008 20101. In the line-by-line radiative 
transfer code used to calculate the cross sections (see above), the broadening pressure 
Pb of the gas species i is calculated from the ideal gas law, i.e. pt, - Ci ■ p instead of 
using a column-density correction to calculate the actual broadening pressure (e.g., see 
|Kasting|[l987| l. However, these differences do not affect the conclusions of previous 
studies. 

The second model verification uses a modern Earth reference profile. Ozone and 
methane concentrations were taken from [Grenfell et aT]p011| l. We run simulations 
with two different thermal radiative transfer schemes. First, we used RRTM ( jMlawerj 
et al.|1997| l, a scheme that is widely used in ID and 3D exoplanet and Earth climate 


studies (e.g., Segura et al.||2003 2005 Roeckner et al.|[2006 Grenfell et al. 2007b 


Kaltenegger et al.|201 1 Stevens et al.|2013l. A second simulation was done with the 


new scheme presented in this work. 

Figure shows the resulting temperature profiles (left panel), the temperature dif¬ 
ference between both models (center panel) and a flux comparison for the upwelling 
thermal flux (right panel). For the latter, we compared calculated model IR fluxes with 
fluxes calculated by the high-resolution line-by-line code MIRART-Squirrl. 

Concerning the temperature profiles, the agreement in the troposphere and the 
lower stratosphere is very good. In the upper stratosphere, small discrepancies can 
be seen, reaching up to about 3 K. However, the overall agreement is good, and such 
small temperature differences are not expected to have a large influence neither on ob¬ 
servables (e.g. spectral signatures) nor on temperature-dependent chemistry. For the 
IR thermal upwelling flux, it is clearly seen that the fluxes do not differ by more than a 
few % throughout the stratosphere. 
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Figure 2: Early Mars: Comparison of temperature profiles from |von Paris et al.H2013b) and this work. 
Carbon dioxide saturation vapor pressure curve as green dashed line. 


Based on these verifications, we consider the new radiative scheme to be applicable 
to further planetary scenarios. 


2.4. Early Mars scenarios 

Simulated model atmospheres were assumed to be composed of CO 2 , H 2 O and O 3 . 
The CO 2 partial pressure is fixed at 1 bar, consistent with upper limits placed on at¬ 
mospheric pressure during the Noachian (e.g., [Phillips et al.|2001[ ). H 2 O is calculated 
according to ambient temperature, assuming a fully saturated atmosphere (i.e., RH=1, 
see [von Paris et al.|2013b| ). This approach is similar to most ID studies of early Mars 
(e.g., Mischna et al.pOOofjColaprete and Toon|2003 Tian et al.|2010| l, but most likely 
over-estimates the amount of atmospheric water. The surface albedo was set to 0.21, 
close to the observed value of present Mars ( jKieffer et al. 11977] !. This value of surface 
albedo allows the model to reproduce present Mars mean surface temperatures. Plan¬ 
etary gravity is 3.73 ms“^, and the orbital distance is set to 1.52 AU. We neglected any 
effect of eccentricity. However, even though Mars’ orbit is eccentric today (e=0.09), 
and eccentricity cycles would occasionally drive the eccentricity to high values (e.g.. 
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Figure 3: Modem Earth model verification: (Left panel) Comparison of temperature profiles calculated with 
RRTM jMlawer et al.||1997^ and MRAC (this work). (Center panel) AT (RRTM-MRAC) between both 
models, (Right panel) IR upwards flux ratio (climate model)/line-by-line. 


Laskar et al.|2004|l, the overall effect would probably be small (at e-Q.2, the mean flux 


increases by only 2 % compared to the circular case). 

The incoming solar irradiation was set to Noachian conditions 3.8 billion years 
ago, i.e. 75 % of today’s irradiation (e.g., Gough 1981| l. The input spectrum was 
taken from Gueymard ( |2004 )l and scaled at all wavelengths accordingly. Note that 
this approach is consistent with previous early Mars studies (e.g., Kasting|1991 Tian 


jet al.|2010] [Wordsworth et al.|2013| ). Other work, primarily on early Earth climate, not 
only scaled the spectrum, but also incorporate the variation of stellar parameters (e.g., 
Goldblatt et al.|2009ll and enhanced UV radiation (e.g., Kunze et al.|2014 1. Enhanced 


UV radiation (see, e.g.. Ribas et al. 20051 of the young Sun could have played an 
especially important role affecting ozone hence the thermal structure of the atmosphere. 

Ozone profiles strongly depend on the assumed oxygen content of the atmosphere, 
the incoming stellar UV radiation, and, in addition, on the catalytic cycles (HO^, NO^:, 
etc.) operating in the atmosphere (e.g., Selsis et al.|2002[ Segura et al.|2007| Grenfell 


etal.|2013| l. Detailed photochemical modeling is warranted for this problem, however 


such relatively complex models would necessitate many more unconstrained bound- 
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ary conditions and parameters (such as surface fluxes of oxidizing and reducing com¬ 
pounds, wet and dry deposition rates, vertical mixing, kinetic data, UV radiation field, 
etc.). Such complex models (e.g., Selsis^etay2002 Segura^e^alj2007| |Hu et al.|2012[ 


Domagal-Goldman et al.|20T4 Tian et al.|2014 Grenfell et al. 2014| l allow for an esti¬ 

mate of possible ranges of ozone concentrations. Motivated by the results from detailed 
photochemical studies, we perform a sensitivity study of the influence of ozone on tem¬ 
perature structure by inserting artificial ozone profiles Co, ■ These profiles are not the 
output of a photochemical model, but are parameterized as a function of pressure p 
based on a Gaussian profile 


Co, ip) = vniTmax ' exp 


(log 


Pmax / 


0.5 


( 2 ) 


where pmax is the pressure of the maximum concentration and vruTmax is the maxi¬ 
mum volume mixing ratio. 

We performed simulations for different combinations of vruTmax and as shown 
in TableIn Fig. we show a sample of the chosen ozone profiles, together with a 
modern Earth mean profile. 


Table 3: Adopted parameter values for eg. [2|in this work. 


Parameter 


Values 


Pmax [mbar] 
vrnrmax [ppm] 


0 . 1 , 1 , 10 
0 . 1 , 1 , 10 , 100 


The ozone layer peak (its location and its magnitude) is the result of an interplay be¬ 
tween two main effects. On the one hand, ozone formation requires photolytic release 
(A < 200 nm) of atomic oxygen (O) e.g. from O 2 , CO 2 , H 2 O, which are usually more 
abundant on the lower levels (followed by fast, three-body formation of ozone). On the 
other hand, UV radiation is more abundant with increasing altitude, hence there is gen¬ 
erally a distinct maximum of ozone concentration at pressures of about 10“^-10“^ bar. 
This motivates the form of our artificial ozone profiles. 

Previous calculations by, e.g., Selsis et ari ( |2002| l and Segura et al. (2007 1 suggest 
a maximum in ozone number density of up to 10'^ cm"^, located at pressures around 
1-10 mbar. They found ozone columns between 10 "^-2 times the modern Earth column 
(Table 3 in |Selsis et al.|2002 and Table 2 in Segura et al. 20071. Both studies focused on 
N 2 - or C02-rich atmospheres, similar to the ones considered here. The oxygen needed 
for ozone formation was either provided by H 2 O or CO 2 . To cover the entire range of 
possible ozone concentrations found by [Selsis et al.| ( [2002| l (i.e., their humid and dry 
CO 2 cases as possible extremes) we varied the maximum ozone concentration between 
0.1 and 100ppm. As shown in Table|^ this leads to relatively thick ozone columns (up 
to 24x the modem Earth values) for a few cases. 
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Figure 4: Sample of chosen ozone profiles. Earth mean profile from [Grenfell et aT]j2Qll) as dash-dotted 
line. 


Table 4: O 3 columns (in units of terrestrial columns, where one terrestrial column corresponds to 315 Dobson 
units) f or the scenarios listed in Table[^ __ 


Pmax [mbar] 

vrnrmax [ppm] 

0.1 

1 

10 

100 

0.1 


0.00023 

0.0023 

0.023 

0.23 

1 


0.00249 

0.0249 

0.249 

2.48 

10 


0.0248 

0.248 

2.47 

24.7 


3 . Results and Discussion 

3.1. Surface temperature 

Figure]^ shows the surface temperature as a function of the parameters in eq. 

For all scenarios, the surface temperatures hrst slightly increase with increasing 
vmr max. due to a slight decrease of planetary albedo and an additional greenhouse effect 
from ozone. For higher vmr max. surface temperatures start to decrease with vmr max (see 
below). 

However, for the Pmax=0.1 and 1 mbar cases, the overall effect is rather small, of 
the order of 1 -2 K. In contrast, for Pmax= 10 mbar, the increase at low virttmax (up to 2 K) 
and the subsequent decrease towards higher values is much more pronounced, with a 
maximum temperature drop of about 8 K compared to the zero-ozone case. 

3.2. Temperature structure 

Figurej^shows the effect of ozone on the calculated temperature-pressure prohles. 
In the middle and upper atmosphere, up to 60 K of warming is seen, depending on 
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Figure 5: Effect of ozone on surface temperature. 


the choice of O 3 parameters. This is of course due to the massive increase in stellar 
heating rates associated with the absorption of solar UV absorption by ozone. Heating 
rates increase by up to a factor of 5 compared to the no-ozone case (not shown). 

Two other effects are readily inferred from Fig. First, for the pinax=10mbar 
cases at vmrmax^lO and 100 ppm, temperature profiles no longer follow the carbon 
dioxide saturation pressure curve (black dashed line in Fig. |^. This suggests that 
carbon dioxide condensation, hence also the formation of carbon dioxide clouds, is 
prohibited throughout the entire atmosphere. For the other cases, the potential cloud 
formation zone is reduced, although the effect is rather minimal at low vruTmax or small 
Pmax- Second, as suggested by the slope of the vmrtnax=100ppm, pmax^lOmbar case, 
the lower atmosphere is no longer convective. This is in contrast to the other scenarios, 
where the lapse rate follows a CO 2 adiabat in the mid atmosphere and a wet H 2 O 
adiabat in the lower atmosphere. 

The strong UV absorption by ozone drastically changes the energy balance in the 
atmosphere. Figure |7] illustrates the effect of ozone on the planetary albedo. The plan¬ 
etary albedo is reduced by up to 40 % at large ozone layers. In terms of additional 
radiative energy deposited into the atmosphere, a reduction of albedo from 0.32 (zero- 
ozone) to 0.2 (maximum ozone effect) corresponds to a forcing of 13 Wm“^. However, 
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Figure 6: Effect of ozone on temperature structure. Reference case (zero ozone) as black plain line. Carbon 
dioxide saturation pressure shown as black dashed line. 


as illustrated in the left panel of Fig. most of this energy is absorbed in the mid¬ 
dle atmosphere. As a consequence of the strong absorption, the incoming solar flux at 
the surface is reduced by about 20%. Therefore, the total radiative flux in the lower 
atmosphere (sum of stellar and thermal fluxes) is reduced, coming closer to radiative 
equilibrium. Hence, less energy is available to drive convection. This then leads to a 
radiative instead of a convective lower atmosphere (as observed in Fig.|^. 

A reduced albedo corresponds to a larger thermal flux that must be emitted by 
the atmosphere to satisfy energy balance and radiative equilibrium at the top of the 
atmosphere. For optically thin atmospheres such as Earth’s, this generally leads to 
an increase in surface temperature. In the case of dense carbon dioxide-dominated 
atmospheres considered here, the thermal emission originates mostly from pressures 
at around 10-100 mbar, since at higher pressures, the atmosphere is optically thick for 
thermal radiation at all wavelengths (see right panel of Fig. [^. For the considered 
scenarios with a substantial increase of local temperature (see Fig. |^, the local IR flux 
also increases. This balances the solar flux deposited in the atmosphere. Most of the 
outgoing thermal flux stems from emission in the 15 //m band of carbon dioxide which 
becomes optically thick at pressures around 10-100 mbar. 
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Figure 7: Effect of ozone on calculated planetary albedo. 


3.3. Discussion 

3.3.1. Oxygen content and catalytic chemistry 

The assumed ozone profiles are dependent on the oxygen content of the early Mars 
atmosphere. G iven that modern Mars ha s only little oxygen in the atmosphere (around 
Yung and deMore||l9991, it seems likely that our high-ozone cases will 


10-3 


vmr, e.g. 


over-estimate the total ozone column by a generous margin. However, the formation 


of ozone also depends on the assumed UV radiation field (e.g., Segura et al.|2007 


Domagal-Goldman a nd Meadows|2()T0 Domagal-Goldman et al. 12014 Grenfell et al. 

2014[|Tian et al.|2014 1 which is relatively poorly known for the young Sun (Ribas et al. 
2005 Observations of younger solar-type stars suggest that far-UV and X-ray emis¬ 
sions strongly decrease with age (e.g.. Ribas et al. 2005|l w hereas the total luminosity is 
generally assumed to increase with age (e.g., Gough|198H Caldeira and Kasting|1992| l. 
Therefore, the ratio between oxygen (hence ozone) production via photolysis (in the 
far-UV) and photolytic ozone destruction (in the near-UV) could presumably change 
dramatically and produce largely varying ozone concentrations (see, e.g., sensitivity 
studies by [Segura et al.|2007 Grenfell et al.|2014 or Tian et al.|2014| l. Therefore, reli¬ 
able quantitative estimates of ozone on early Mars are challenging to obtain. Detailed 
photochemical models would probably only allow for an order-of-magnitude estimate 
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Figure 8: Effect of ozone on downwards solar (left) and upwelling thermal fluxes (right). Captions coiTe- 
spond to both panels. 


of such quantities, since boundary conditions (such as wet deposition of hydrogenated 
species and oxidation of reduced volcanic gases) are essentially unconstrained. 

In addition, ozone chemistry heavily depends on the operating catalytic cycles in 
the atmosphere (e.g., Crutzen|1970 Grenfell et al.|2013] l. On Earth, most of the cycles 
are related to NO^: and HO^:. Possible NO^ sources on early Mars would be cosmic 
rays (e.g., Jackman et al.|1980) l or lightning (e.g., Segura and Navatro-Gonzalez 2005| l, 
but these are hard to estimate. It is unclear whether significant HOx cycles could op¬ 
erate on early Mars, since calculated stratospheres are very dry (see Fig. |^, with vmr 
about 10“^-10 This is up to several hundred times drier than in the modern Earth 
stratosphere (H 2 O vmr of a few times 10 ® up to 10“^). Note, however, that HOx cycles 


are responsible for the low oxygen content in the present Mars atmosphere (e.g., Nair 


|et al.|[T994l [Yung and deMore]|1999[ [Stock et^|2012a]bl l since they quickly recycle 
the products of CO 2 photolysis (COh-O) back to CO 2 . Without such recycling, the 
martian atmosphere would be more 02-rich than today, probably at around 5-10% O 2 . 
CO 2 would still remain the dominant atmospheric species (e.g., Nair et al.|1994| l. This 
can be explained by the fact that at increasingly higher O 2 concentrations, O 2 UV ab¬ 
sorption starts to shield the CO 2 from becoming photolyzed since O 2 photodissociation 
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cross sections are much larger than corresponding CO 2 cross sections (see, e.g., Yung 


and deMore||1999 Selsis et ar]|2002|l. Even on Venus with a very dry stratosphere, 


HO^ cycles operate, using HCl photolysis as source for H ( |Yung and deMor^|1999| l. 
As an alternative destruction cycle, SO;^ cycles could be possible, since higher amounts 


of sulphur-bearing species are expected in the early Mars atmosphere (e.g., Farquhar 
et al.|2000 Halevy et al.|2007 Johnson et al.|2008 [Ramirez et al.|2014) l. 


Investigating these possibilities is however beyond the scope of this sensitivity 
study that focuses on climatic effects. 

3.3.2. Eddy diffusion and convection 

As a result of the suppression of convection in the lower atmosphere (see Fig. for 
the high-ozone case, the atmosphere becomes stably stratified and vertical transport is 
significantly reduced. This has a potentially large impact on (photo)chemistry in such 
atmospheres. In general, photochemical networks include photochemical reactions (in¬ 
cluding e.g. 3-body reactions, photolysis, decomposition reactions, etc.) as well as a 
vertical mixing, which is parameterized by so-called Eddy diffusion. For the convective 
troposphere, ID photochemical models of Mars or Earth use constant, relatively high, 
Eddy diffusion coefficients (e.g., |Massie and Hunten|19^ Yung and deMore|1999j l. In 
a stably stratihed atmosphere, the diffusion is most likely less efficient. The magnitude 
of the effect on convection and transport depends, among others, on the topography. 
Thus, full 3D dynamic simulations are probably needed to investigate the influence on 
chemistry. This is, for modern Earth and Earth-like (terrestrial) exoplanets, mostly 
important for biogenic trace gases such as methane and nitrous oxide. In the lower at¬ 
mosphere, these are more strongly influenced by transport rather than chemistry (e.g., 
Segura et al. [2003] [2005[Grenfell et al.|2007a|b| ). 


Exploring this issue, while beyond the scope of this work, would however be inter¬ 
esting in the context of fully coupled climate-chemistry simulations. 

3.3.3. Water loss from early Mars 

The very thin atmosphere of current Mars is likely shaped by significant atmo¬ 
spheric loss. Isotopic signatures of atmospheric escape have been found for nitrogen 
and oxygen (e.g., Jakosky and PhiIIips|200I Fox and Hac|20I0 Gillmann et al. 201 Ij l. 
Recent D/H isotopic ratio measurements of the Curiosity rover on Mars suggest that 
hydrogen escape and water loss has been ongoing over the last few billion years (e.g., 
[Mahaffy et al.|20f5] l after the (presumed to be dense) Noachian atmosphere has been 
eroded. Modeling the escape of water is a very important issue for atmospheric evo¬ 
lution and habitability (e.g., Kasting|ri988 Kulikov et al.||200^ . The water loss rate 
depends on the stratospheric water content, where water is photolyzed and subsequent 
loss of H atoms occurs (e.g., Kasting|1988 [Kasting et al.|1993] l. 

Figure shows the calculated H 2 O prohles. Note that we assume a relative hu¬ 
midity of unity throughout the atmosphere hence H 2 O concentrations are likely to be 
over-estimated in the mid- to upper atmosphere. On Earth, the stratospheric H 2 O con¬ 
tent is mainly controlled by the efficient cold trap near the tropopause (around 10 km 
altitude). With increasing temperatures due to the ozone layer, water is efficiently 
trapped in the troposphere. This leads to the dry stratosphere and, consequently, low 
water escape rates. 
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Figure 9: Effect of ozone on calculated water profiles. 


For the early Mars scenarios considered here, cold traps develop even in the absence 
of temperature inversions. However, given that surface temperatures are very low, and 
the cold trap is located high in the atmosphere, stratospheric H 2 O concentrations are 
very low. This suggests that water loss might be very slow on early Mars. However, 
when increasing the ozone concentrations in the model atmospheres, the altitude of the 
cold trap changes dramatically, due to the change in the temperature structure (see Fig. 
1^. This leads to a significant increase of stratospheric H 2 O content (up to two orders 
of magnitude), and possibly impacts the water escape rates. Hence, results suggest that 
water escape rates might be sensitively influenced by the ozone content. Note, how¬ 
ever, that an increase in stratospheric H 2 O concentrations could also lead to enhanced 
destruction of ozone through catalytic HOjc chemistry (see, e.g., the parametric simu¬ 
lations of |Tian et al.|2014) l. Therefore, changes in cold-trap locations might limit the 
increase of ozone, hence limit also H 2 O escape. To further investigate this issue would 
need consistent climate-chemistry models to assess this feedback between vertical tem¬ 
perature structure, stratospheric water content and ozone chemistry. 

To estimate the amount of water lost via atmospheric escape, we adopted the same 
approach as Wordsworth and Pierrehumbert ( 20T4) l. Escape was assumed to be diffusion- 
limited and controlled by the cold-trap temperature T^t and cold-trap water concentra¬ 
tions. The diffusion coefficient of H 2 O through a CO 2 atmosphere was taken from 
approximative equations of Marrero and Mason (19721. 

These approximations lead to water escape fluxes of, even under the most favorable 
escape conditions, only about 1 cm of water lost within 10® years (1 Gyr). This has to be 
compared to the actual atmospheric water column of about a few 10 m of water in our 
simulations. Furthermore, the original Martian water inventory is estimated to be much 
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larger, about 10^-10^ m (e.g., McKay and Stoker 1989 Lammer et al. 2013 Lasue et al. 


2013| l, Currently, most of the water is thought to either reside in the cryosphere or to be 
bound to minerals in the crust (see discussion in, e.g., [Lammer et al.|2013|and|Lasue 
jet al.|2013» . 

In addition, we point out that current best estimates of the hydrogen escape flux 
suggest about 1-2 m of water lost in 1 Gyr (e.g., Chassefiere and Leblanc|2004 Zahnle 
2015| l, about 2 orders of magnitude more than we calculate. This is due to the fact that 
molecular and atomic hydrogen are the main hydrogen-bearing species in the upper 
Martian atmosphere, not H 2 O, as we assumed above. H and H 2 have concentrations 
about 2 orders of magnitude higher than H 2 O (e.g., Nair et al.|1994 Yung and deMore 
1999| l. H and H 2 are derived from H 2 O photolysis, and since water is controlled by 
surface conditions and constantly replenished, H and H 2 can build up until concen¬ 
trations reach an equilibrium between the photolytic source and atmospheric escape. 
Addressing the connection between cold-trap location, thermal structure and H 2 O es¬ 
cape necessitates coupled climate-chemistry simulations, which is a potential direction 
of future research. 


3.3.4. Spectral appearance 

In an exoplanetary context, an interesting question with respect to ozone is the 
possibility of detecting it remotely with spectroscopic observations. It has been pro¬ 
posed as a so-called biosignature gas (e.g., Leger et ar]|1996 Selsis||2000 Schindler 


and Kasting|2000 Des Marais et al.|2002| Selsis et al.|2002 1, and much work has been 


invested into atmospheric and spectral modeling regarding detection and interpretation 


of ozone signatures (e.g., Selsis et al.|2002 Segura et al. 2003 2007 Kaltenegger et al. 

2007 

Kaltenegger and Traub 2009 

Rauer et al.|201 1 Hedelt et al.|20 13 

Grenfell et al. 

2014 

1. 


Considering early Mars as a potential exoplanet. Fig. 10 shows the resulting emis¬ 
sion spectra for five selected cases. They have been calculated with the same high- 
resolution line-by-line code used for the model verification ( Schreier and B6ttger|2003| l 
and were then binned to a coarse spectral resolution of R= ^ = 100. As pointed out by 
Selsis et al. ( 2002| l, the ozone band at 9.6 pm coincides with a carbon dioxide hot band 
and lies in the wings of the strong dimer absorption feature (see, e.g., [Wordsworth et al. 
2010 for the most recent carbon dioxide continuum absorption data). Therefore, except 


for the strongest ozone column, the shape of the spectral band is not greatly sensitive 
to atmospheric ozone amount. 


To illustrate this. Fig. 11 zooms in on the spectral region around the 9.6pm fun¬ 
damental band of ozone. For comparison, we also show line-by-line calculations with 
removed ozone (marked ”no 03” in Fig. 0 to investigate the effect of ozone on the 
spectral appearance and a possible masking by CO 2 . It is clear that an effect is only 
seen for the highest ozone columns (higher than the modern-day terrestrial column. Ta¬ 
ble!^. To detect this, relatively high spectral resolution (here, R=100) and reasonable 
signal-to-noise ratios (of the order of 4) are needed. These values are not achievable 
with currently planned instrumentation even for Earth-sized or larger planets, let alone 
a planet the size of Mars (see, e.g., Hedelt et al.|2013 1. 

The strong carbon dioxide fundamental band around 15 pm is very sensitive to 
middle-atmosphere temperatures (see Fig. [^, and the effect of the changing temper- 
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Figure 10: Emission spectra (spectral resolution R=100) for selected cases. 


ature structure is clearly apparent in this band. In the hypothetical case that all stel¬ 
lar parameters (including UV radiation field) were known, and consistent atmospheric 
modeling were used, then this band would constimte an excellent indirect probe of at¬ 
mospheric ozone content. We note, however, that other radiative gases could also be 
responsible for a temperature inversion (e.g., SO 2 , CH 4 , etc.). This then is an important 
challenge for retrieval algorithms (e.g., |von Paris et al.|2013c [Irwin et al.|2014| l. 

In addition to that, the broad 6.3 jum band of water, superimposed by part of a 
carbon dioxide dimer absorption feature, is also sensitive to middle-atmosphere tem¬ 
peratures and thus would be an indirect probe of ozone content. 

Summarizing, in the cases considered for this work, ozone is hardly seen directly 
in the spectra, except at high ozone contents. Indirectly, however, due to its impact on 
the thermal structure, the effect of ozone is clearly seen in the spectra. 

As outlined in jKitzmann et al.jpOl Ij l or jvon Paris et al.| ( |2013d| l, we use the model 
output of the stellar radiative transfer code to produce low-resolution albedo spectra. 
In Fig. 12 we show the UV-visible part of the spectram. For the low-ozone reference 
case, the spectral albedo clearly follows the Rayleigh scattering slope expected for 
dense carbon dioxide atmospheres. The presence of ozone leads to a strong reduction 
of albedo in the Hartley (0.24//m< /l<0.34/im) and Chappuis (0.4//m< A <0.7jum) 
bands. Additionally, the location of the ozone maximum also influences the spectral 
albedo in the Hartley band, as shown in Fig. 


12 For cases where the maximum is 


located at higher pressures, the UV radiation penetrates deeper into the atmosphere 
and thus is more susceptible to Rayleigh scattering by carbon dioxide. 
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Figure 11: Zoom in the 9.6/tm band of ozone. Effect of including ozone in spectral calculations (R=100). 


4. Conclusions 

Ozone is a radiatively important atmospheric species due to strong UV absorption 
bands. It can significantly change the atmospheric temperature structure and thus af¬ 
fect, e.g., cloud formation or atmospheric transport. 

We have investigated the influence of ozone on the atmospheric temperature struc¬ 
ture of early Mars. We simulated a 1 bar carbon dioxide atmosphere with different, 
fixed ozone concentration profiles. To do this, we have developed a new IR radia¬ 
tive transfer scheme that incorporates molecular absorption of carbon dioxide, water, 
methane and ozone. Calculated IR cross sections cover a wide range of temperatures, 
pressures and concentrations. Hence, this scheme can be applied to a wide range of 
possible planetary scenarios. 

Results suggest that the impact on surface temperatures at small to moderate ozone 
concentrations is probably not large, of the order of 1-3 K, depending on the overall 
ozone column. The increase in surface temperature is mostly due to the decreased 
albedo, since ozone strongly absorbs the incoming UV and visible solar radiation. For 
high ozone concentrations, surface temperatures drop by up to 8 K due to a change in 
energy balance. In the upper and middle atmosphere, temperatures increased by up 
to 60 K upon introducing ozone. The resulting increase in thermal flux balances the 


19 



















wavelength [/um] 


Figure 12: Effect of ozone on spectral albedo. 


radiative forcing of UV absorption and leads to the observed surface cooling at high 
ozone concentrations. 

As a consequence of the stratospheric warming, the cloud forming region is re¬ 
duced. In the case of a thick ozone layer (comparable to or larger than the terrestrial 
ozone layer), cloud formation is inhibited completely. 

For large ozone columns, convection is inhibited. Instead model atmospheres are 
fully radiative due to a strong reduction of incoming solar flux in the lower atmosphere. 
This probably decreases the vertical transport. 

Upon increasing the ozone content in the model atmospheres, stratospheric water 
content strongly increased, due to a change in thermal structure and a change in cold 
trap location. Compared to zero-ozone simulations, water concentrations increased by 
up to two orders of magnitudes. This is important for assessing possible atmospheric 
water loss on early Mars, but could be a self-limiting effect due to enhanced catalytic 
HO^ cycles that destroy ozone. 

Future work aims at performing fully coupled climate-chemistry simulations of 
early Mars as well as carbon dioxide-dominated atmospheres near the outer boundary 
of the habitable zone. 
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